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Abstract. 

We propose alternative determinant representations of certain form factors and 
scalar products of states in rational Gaudin models realized in terms of compact spins. 
We use alternative pseudo-vacuums to write overlaps in terms of partition functions 
with domain wall boundary conditions. Contrarily to Slavnovs determinant formulas, 
this construction does not require that any of the involved states be solutions to the 
C$ • Bcthc equations; a fact that could prove useful in certain non-equilibrium problems. 

Moreover, by using an atypical determinant representation of the partition functions, 
we propose expressions for the local spin raising and lowering operators form factors 
which only depend on the eigenvalues of the conserved charges. These eigenvalues 
^ 1 define eigenstates via solutions of a system of quadratic equations instead of the usual 

C^l . Bethe equations. Consequently, the current work allows important simplifications to 

| numerical procedures addressing decoherence in Gaudin models. 

o 

1. Introduction 

Integrable models based on the generalized Gaudin algebra [TJ [2] have, in recent years, 
found a large ensemble of physical applications ranging from the mesoscopic BCS model 
[31 HI E] to the central spin Hamiltonian ^ [H |gl |S1 HOI IHl H21 H31 HH HSl HSl 113 US] 
through a variety of cavity based constructions relevant for quantum computing 
proposals [191 EQl [21]. The fact that their integrability does not necessitate strong 
restrictions on the model's parameters also makes them a remarkable playground to 
study externally tunable physical systems. 

The exact eigenstates of Gaudin models are obtainable by finding sets of complex 
parameters (rapidities) which are solutions to an ensemble of non-linear algebraic 
equations known collectively as Bethe equations. However, the efforts to numerically 
solve these equations in a systematic fashion have shown it to be a challenging task 
|22[ 1231 EU [25] . Recently an important improvement [261 [27J [28] has been achieved by 
exploiting a non-trivial change of variables based on the correspondence between Bethe 
equations and ordinary differential equations [29, 30J. In doing so, one can rewrite the 
problem in terms of quadratic equations depending on a new set of variables A(e$) which 
are directly related to the eigenvalues of the model's conserved charges. 
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Using Lagrange's polynomial basis it was possible to implement an approach 
allowing one to extract the rapidities from a given set of A(e$) [27J. In doing so, 
one could turn to Slavnov's determinant [31] in order to compute scalar products 
and local operator form factors which are the elementary building blocks needed to 
address physical quantities. However, this work also motivated the search for simple 
representations of these quantities expressed, not in terms of the rapidities themselves, 
but directly in terms of the easily found A(e;) variables. The current paper addresses 
this question and proposes to do so by using a non-standard determinant expression for 
the partition function with domain wall boundary conditions. In conjunction with the 
existence of two distinct representations for the eigenstates we find single determinant 
expressions for overlaps and spin raising/lowering operators form factors. 

The paper is organized as follows. In section |2] we briefly review the Algebraic 
Bethe Ansatz (ABA) as applied to Gaudin models. Putting the emphasis on the two 
possible quantization axis ±5, we find a simple transformation between two equivalent 
representations of any eigenstate of the system. In Section [3] we then derive a 
determinant expression for the partition function with domain wall boundary condition 
which is used in Section H]to write scalar products of Bethe states as simple determinants. 
Section [5] concentrates on deriving determinant expressions for the form factors of local 
spin operators. In Section [6] we discuss two possible applications of the obtained results 
to non-equilibrium problems. 

2. Algebraic Bethe Ansatz 

Let us first introduce the generalized Gaudin algebra defined by the operators 
S x (u), S y (u), S z (u) satisfying the commutation relations [TJ [2]: 



[S x (u), S»] = i(Y(u, v)S z {u) - X(u, v)S z (v)), 

[S y {u),S z {v)] = i(Z(u,v)S x (u) -Y(u, v)S x {v)), 

[S z {u), S x {v)} = i{X{u, v)S y (u) - Z(u, v)S y (v)), 

[S»,S»] =0, K = x,y,z, (1) 

where it, v 6 C. In this paper, we will deal only with the rational family of Gaudin 
models for which 

X(u, v) = Y(u, v) = Z(u, v) = (2) 

it — v 

For a given number of excitations M, the ABA allows one to find eigenstates of the 
transfer matrix T{u) = S 2 (it) using the following construction 

M 

lAi.-Aj^) =US + (Ai) |0). (3) 

i=l 
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Here S + (u) = S x (u) + iS y (u) are generalized creation operators parametrized by the 
complex variable u. The pseudovacuum |0) is defined as the lowest weight vector, i.e. 
S"(«) |0) = 0, V w G c. 

States of the form become eigenstates of 

T(u) = S 2 (u) = X - (S + (u)S-(u) + S"(w)S + (u) + 2S»S») (4) 

provided the M rapidities Aj are solution of a set of coupled non-linear algebraic 
equation: the Bethe equations. For rational models, these equations can be written, 
in general, as 

M 



A, — A, 



with 



S*(A«) |0) = F(Ai) |0) (6) 

defining the lowest weight function F(u). 

Since one can show that [S 2 (w), S 2 (t> )] = 0, the operator-valued residues {Ri...Rn} 
of S 2 (it) at its arbitrarily chosen poles u G {ei, ...ejv} 

allows one to define a set of A" commuting hermitian operators R{. These 
become constants of motion for any integrable Hamiltonian obtained through linear 
combinations using coefficients rji G R: 



N 



H = Y J mRi- (7) 



i=i 



2.1. Correspondence between pseudo-vacuums 

When dealing with Gaudin models realized in terms of operators bounded from above 
and below, we have the freedom of defining the ABA using either the ±5 quantization 
axis. Including an external magnetic field -z, the two realizations in terms of A" local 
su{2) spin operators of lenght \Si\ are given by: 



|o) = ||...|) |0> = |t-t> 

N N 

i=l " tl i=i U t% 

u — e» ^— ' u — €i 

i=l i=l 

AT „_ iV 

S» = £^^ = C(u) S-(n) = E ^ 



U — €i — ' U — Ei 

i=l i=l 
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where j" (1-) respectively represent the highest (lowest) weight state for each local spin. 
Note in passing that this readily excludes any model containing bosonic degrees of 
freedom such as Jayne-Cummings-Dicke-like models. Nonetheless, for any realization in 
terms of finite magnitude spins or pseudo-spins, both constructions are available. 
The generic states containing M up spins 

M 

\\ X ...X M ) = J]B(Ai)|4-...|> 



iml-.^-m) = n it - 1) (9) 

i=i 

turn into eigenstates of the transfer matrix provided the rapidities Aj or \a satisfy the 
Bethe equations ([5]): 



AT i „ , „ M 



^) = -E^4= E ^ 

k=i K * v i=i(^0 J 



TV | „ , „ N—M 



= -E—- 1 = E > ( 10 ) 

while the eigenvalues of S 2 (u) are then given by 



^)] 2 -^»- 2 E^ + E^t E^ 



^) - - >h - 2 E ^ + E ^ ( E z± n 

:n) 



— ' U — Hi £ — ' U — Uj \ *■ — ' U — IX i 

=1 I* i=l ^ \i=l(^i) PJ 



The poles of these eigenvalues at u — €j give the eigenvalues of the commuting 
operators Ri, which are themselves read off from the poles of the S 2 (u) operator. 
Specializing to the non-degenerate case (e^ ^ tj V i ^ j), we find: 



2Sf ^ 2&-& r A A 2 2 A 21 5, 



«. J = -^+ E - Fh = -Er4r + 7+ E 



i=l(=#) J 1 11 i=l J y i=l(=#) 3 

2Sf ^ rf 2 2 A 215. 



9 i^ti-tj \Si\ jri^i-N 9 ^ 



J I 



(12) 

Unsurprisingly, one has the same conserved charges R} = Rf. In order to find 
a transformation leading from one representation of a given eigenstate to its other 
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representation, it is sufficient to insure that every eigenvalues rj are the same in both 
cases. In doing so, one easily sees that the transformation 

A"(eO = A A ( ei ) - - (13) 
9 

does give the correspondence between both representations of a given eigenstate. Here 
we introduced the variables 



A 



A/ \ \ ^ 



M 1 



N—M 



= E 7^7' ( 14 ) 

which are directly related to the eigenvalues of the commuting Gaudin Hamiltonians 
Ri (see (TT2|)). 

One should keep in mind that the transformation is exclusively valid for states 
which are solutions to the Bethe equations (eigenstates) and that, evidently, the two 
representations can still differ by a normalization factor. Moreover, one should note that 
A(e«) are sufficient to allow a direct construction of the eigenenergies of any integrable 
Hamiltonian of the form H = YliLi Vi^i with r)i G K. 

Working with the rapidities {\x-.-\m}, {^x---^n-m}i establishing a transformation 
between both representations would only be possible by solving a further set of non- 
linear equations whereas here, using the A(e,)'s, it is remarkably simple. 

2.2. Bethe equations for A(q) 

As briefly mentioned in the introduction, the A(ej) variables provide an extremely useful 
representation of the eigenstates in the sense that they obey a set of algebraic equations 
which is much simpler than the underlying Bethe equations obeyed by the rapidities \. 

For simplicity, the remainder of this paper will focus on non-degenerate realizations 
in terms of spin \ operators = |). It was shown [32] and exploited numerically 
|26[ 127] that, in this case, solutions to the system of N quadratic equations: 

[A „ fe)]2 = ^ A-fe) - A-fc) _ 2 (1B) 

i=l(^j) 3 6t 9 

are in one to one correspondence to solutions of the Bethe equations fflU]) via the 
definitions (fl4l) . It is a trivial matter to verify that transformation (|T3|) is consistent 
with both versions of eq. (TIB]) . 
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3. Partition function 

Due to the relative simplicity of solving eqs (fl5l) . it becomes highly desirable to be able 
to access physical quantities in terms of simple expressions involving exclusively the 
A(ej) variables. While Slavnov determinants fulfill such a role in terms of the rapidities 
Aj, in the rest of this paper we will derive determinant expressions for scalar products 
and form factors of local spin operators in terms of the A(ej) variables. 

The first step, carried out in this section, is to show that the overlap of a generic 
Bethe-like state ([9]) with an "infinite magnetic field (g = 0)" eigenstate (le^.-.e^) = 

M 

Y\ |4 ••• 4)) is writable as: 

3=1 

{e ix ...e iM |Ai...A M ) = DetJ 

M l 

> 2 — — A ^ a = b 

1,1, - \ c=1( _, a) e ic (lb) 



In order to show this, one can start from the explicit construction of the state 
\\i...X M ) (eq. ©), which leads to the formal expression: 



M 



Here {P} is the ensemble of possible permutations of the indices {ii-.-Uf} and Pi denotes 
the i th element of the given permutation. By isolating in ( TT7|) the terms which depend 
on Am, one finds that the overlaps obey the simple recursion relation 

M 1 

( e U--- e «A/ l^l—^Af) = / , T {tii—hj—tiM |Ai...Aa/_i) , (18) 

3=1 3 

where ej 1 ...ej i ...ej M ) is the state with M —I excitations, for which has been removed 
from the ensemble {e^.-.e^}. 

This is obviously a rational function of Am, which goes to zero when Am — > oo 
and has only simple poles at every Am = To show that it does obey the recursion 
relation, it is therefore sufficient to show that the proposed determinant representation 
(fTEj) has the same poles Am = ^ and the same residues (e^.-.e^ ...€i M |Ai...Am-i) at 
these poles. 

The determinant in (fTEl) clearly only has single poles at Am = Q r Indeed, the 6j. 
pole comes only from the diagonal element Jjj which, via — A(e^), contains the term 
Xa{ 1 _ e . ■ The residue is trivially given by the determinant of the minor obtained by 
removing line and column j after taking its Am — > Q limit: 
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lim (\ M -e ij )'DetJ = 'DetJ j (19) 

with 

{M 1 M-1 l 1 
^ ^ , a = b (a,b ^ j) 
c= j^ a) ^-^ ti e ^- Xk 6i «- €i i (-20) 
a^b(a,b^j) 

M x 

The diagonal elements of this matrix evidently reduce to 



M-1 .| 

> — and therefore correspond to the representation (fl6|) of (e^...^. IAi...Am-i) 

^— ' ' , — A ' 



-la xv « 



proving the determinant obeys the recursion relation ([18]) . 

Verifying that, for a single rapidity Ai, the projection (e^ |Ai) = : — is indeed 
equivalent to the 1 by 1 version of the above determinant (— = — e , _ A ) then 
completes the proof. 

This construction is in fact nothing but the partition function with domain wall 

boundary conditions which one would obtain using a reduced model which contains only 

s+ 

the M states excited in the left state, i.e. using operators 5(A) = J2j=i '■ 



M 



[e il ...e iM |Ai...Am) — (tiitia ■•• Ti^l XX-^(-Ai) |>Ui>Ua ■•• U M ) ■ ( 21 ) 



i=l 



Expression ( fl6l) can however be contrasted with the appropriate limit of the more 
frequently encountered Izergin (33J, [3U [35] determinant representation of such a scalar 
product, i.e.: 



M 



IAi-Am) = -r 4 ^ T , DetK 



M M 

n - a ) n < 

i>j=l j<k=l 



K ab = prj. (22) 

K - Aa) 2 

which is not simply writable in terms of A(e$). One should keep in mind that the 
determinant expression (fl6|) (just as ff22l) ) is valid for any set of complex parameters \ 
and does not require them to be solution to the Bethe equations. 

Finally, it is worth pointing out that due to the invariance under the exchange of 
the sets {e^.-.e^} and {Ai...Aa/} (as evidenced by expansion ffl7]) ). one could also write 
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the projection in terms of the rapidities themselves as the determinant of the following 
alternative M by M matrix: 

M ^ M 1 

■U={ J^ a) *a-K tiK~e ic . (23) 

4. Scalar products 

The scalar product between two generic states (eq. Mj built out of the two different 
representations using respectively M and N — M rapidities is then writable as 

N—M M 
i=l j=l 

= (t ... t h-"AT> , (24) 

where {ui, ...vn} = {^[■■■^'n-m} U {Ai-.-Ajv/} is the union of both sets of rapidities and 
has cardinality N. In doing so, we are once again dealing with a partition function with 
domain wall boundary conditions, this time using the full set of N local spins. The 
results of the previous section are directly usable and lead to the determinant of the 
N x N matrix: 



(fi[...fi' N ^ M |A x ...Am) = ~DetK 



K„h = < , . e a — e, 

c=l(^a) 



N 1 



a ^ b 



f-a—f-b 
N 1 

£ k\e a )-W\e a 

^— ' e„ — e r 

C =l(^a) 



(25) 

We note that for any ensemble of rapidities whose union has cardinality ^ N, both 
states would have different magnetizations and would therefore be orthogonal. 

Contrarily to the traditional Slavnov determinant for (A^.-.A^ \\i---\m) which is 
only valid when one of the two states is a solution to the Bethe equations, the current 
expression has no restriction on any of the two sets of rapidities. Provided the //-state 
is an eigenstate, it corresponds to an alternative A'-state using transformation ffT3]) and, 
in this specific case, we have 
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{\' 1 ...\' M |A x ...Am) oc (^'l—^jv-Af |Ai...A M ) = DetK 

- 1 2 
^ = S c= 7^ a) e - " £c • ( 26 ) 



i 



While the issue of the normalization will be discussed in the next section, we showed 
that by mixing both representations one can write the scalar products of unnormalized 
states in terms of A(ej) variables. 

4-1. Normalization 

For any state which allows both representations |Ai...Am) or |jUi.../ijv_jy), the actual 
norm of either representation expressed in terms of the A(ej) variables remains elusive. 
However, their scalar product (/ii.../ijv_M \\\...Xm) is straightforwardly writable as a 
determinant. Since both representations correspond to the same normalized state 
\^i—^M) Norm = ■fi-\ti>i~>fJ'N-M) = jfc |Ai...Am), the mixed representation allows us 
to write 

N-M M 

N,N X = (| ... tl II Bi^JlBiXi) \i ... I) = DetG (27) 

1=1 8=1 

with the N by N matrix given by 



t4 A A (e„)-A^e„ 

r ab = \ c=l(^a) 



G nft = < ,=7^ ea_ec • (28) 



In the specific case of eigenstates of the system, the correspondence f[T3"j) allows us 
to write it as 

f E "2A A (e a ) + - (a = 6) 

Provided expressions for the form factors (/x' 1 ... / u^ v _ M | (9 |Ai...Am), this product is 
sufficient to write the eigenbasis representation the O operator: 



f= E 

{x' 1 ...x' M },{x 1 ...x M } 





■■^'n-m\ 


o\ 


\\...\m) (i*>i~-Hn-m\ 


{Hi... fJ>N-M 


\\...\m) (fA—f^N-M 


X' 1 ...X' M ) 



(30) 
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Here, one should understand that the notation uses the following correspondence 
^\fi 1 ...fjL N - M ) = ]^|Ai...Am) and ^ \fi[...fi' N _ M ) = j^j \X[...\' M ) while the double 
sum covers twice a full set of eigenstates. 

For any state, be it an eigenstate or not, which is writable using both 
representations, expectation values of a given operator would also be normalizable by 
writing them as: 



(fJ>l---(J>N-M\ 




(/iX.../i/V-M 


Ai...Am) 



(31) 



Having even shown how to go from one to the other via the transformation f|T3|) . 
we know with certainty that both representations are available for eigenstates of the 
system. However, for a generic state built out of arbitrary rapidities {Ai...Aj\/} it is not 
assuredly possible to build an equivalent {//i.../^m} representation. Still, in Section l6\2l 
we discuss a possible scenario where, without being an eigenstate of any given static 
model, a physically relevant time-dependent state would be such that these two possible 
representations exist at any time making ( 13TI) a usable construction. 

5. Form factors 

In this section we derive determinant representations for form factors of local spin 
operators. 

5. 1 . Sf form factors 

The solution to the quantum inverse problem for the models considered here allows one 
to write local spin operators in a remarkably simple fashion. Indeed, local spin raising 
operators are simply given by: 



Sf= lim(7-€i)S( 7 ). (32) 

This fact allows one to derive simple expressions for their form factors. Using the 
multi-representation construction, we obtain for the form factor between unnormalized 
states with M and M+l up-spins: 



(^[■■■^'n-m-i] &t \X\---Xm) — (Ai...A M | S i |/4---Mjv-M-i) 




= lim(7 — €i) detJ, (33) 

where the matrix J is given by eq ( |25i) with the values of A I/ (e a ) obtained for the 

ensemble {/i' 1 .../i^ v _ M _ 1 , 7, Ai...Am}- The determinant has a single pole at 7 = and 
consequently, since 



On Gaudin model's scalar products and form factors 



11 



lim A'fo,*) = A'fo) + A A (e,) + 

€j — €i 

lim( 7 -e J )A^(e i ) = -l, (34) 

the resulting form factor is simply given by the determinant of the (N — 1) x (A — 1) 
matrix: 

(^'i-A*jv-Af-i| S t IAi-.-Aat) = detJ' 

AT 



J2 — A^'(e a )-A A (e a ) a = b & z) 



J^^~^ ■ ( 35 ) 



i 



a ^ 6 i) 



which excludes from the sums as well as line and column i. 
5.2. Sf form factors 

The Sf form factors are obtainable in a similar fashion except for the fact that one 
needs to explicitly use commutation relations of A(u) and B(u) operators. The inverse 
problem gives us 



5? = -lim (T-eM7), (36) 



and therefore 



N-M M 

lim( 7 - Q )(t,...,tl II BQQAWUbMII,...,)). 



1=1 j=l 



(37) 



Using the commutation relations (jTj), it is a straightforward exercise to commute 
the A operator until it reaches the right and acts on the pseudo-vacuum \l, In 
doing so, one obtains the following sum: 

(v'v-v'n-mI ^1 |Ai...A M ) 

i , A i , 



fi 1 ...fi N _ M 



Xi...\ M ) + V] r (^'i-a4-m| s t Xi—Xj—Xm) ■ 

— U — An I 

(38) 



where every term is writable as a determinant. However, we did not manage to reduce 
this sum to a single determinant. Such a feat is possible [36] for (X! X ...X' M \ Sf \X\...Xm) 
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using the Slavnov construction in terms of the rapidities since all determinants then 
differ by a single column. Consequently, it appears that the particular expression found 
here cannot be useful in any numerical application which involves the computation of 
a large number of S z form factors; even more so considering the fact that it would still 
require explicit knowledge of the rapidities Xj. While obtaining rapidities from the set 
of A(ej) is possible following the procedure outlined in [27], having done so would clearly 
make the use a single Slavnov determinant a better suited approach to the computation 
of the form factors. 

Nonetheless, this construction still has the advantage that, contrarily to Slavnov's 
formulas, it remains valid even when both {//} and {A} are not solutions to Bethe 
equations. In Section I6.2[ we discuss a potential scenario in which one could explicitly 
exploit this fact. 

6. Applications 

6.1. N on- equilibrium dynamics 

One of the central motivations behind this work was to numerically address the 
decoherence in the central spin model. It describes a central spin So coupled to an 
external magnetic field Bz and interacting via non-uniform hyperfine couplings A, with 
a bath of N spins Sj. Its Hamiltonian is obtained using a single integral of motion 
H = ^Rq and using the correspondence B = —-, eo = Aj = — — which leads to: 



In order to compute the non-equilibrium dynamics of a generic initial state writable 
as Bethe-like construction one can use the set of determinants proposed in this work 
and alleviate the necessity of explicitly finding rapidities Aj in order to describe the 
eigenstates. Starting from an initial condition given by a coherent superposition of the 
central spin and any arrangement of the bath spins with the spins {ii-.-im} pointing up 
and the rest pointing down: 



N 




(39) 



i=i 



1^(0)) = a |to; I ... tii - ti M - I) + a \U; I - t<i - ti M - I) 

= oc |e ; e it ...e iM ) + f3 \e ix ...ei u ) , 



(40) 



one can write the coherence factor as: 



(mis, 




\{\}„) ({fi}n\ S l"Mm) {{lAm h 




n,m 



(Wn \{Mn) (Mm |{A} m ) 



(41) 
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where m, n respectively cover the full sets of M and M + 1 excitations eigenstates 
with energies o; m n . In light of the work presented here it should be clear that the 
eigenenergies, the form factors and the overlaps of the initial condition with eigenstates 
are all writable exclusively in terms of A(ej) variables. The proposed expressions become 
particularly useful for intermediate system sizes such that the extra computational cost 
associated with N by N determinants (instead of M by M for Slavnov's formulas) 
outweighs the cost of extracting the rapidities A from the set of A(e;). 



6.2. Dynamical Bethe Ansatz 

Finally, considering that a dynamical Ansatz |Ai(i)...Ajw(t)) can > i n certain scenarios, 
describe exactly the non-equilibrium wavefunction for Gaudin models [37J , the ideas 
developed in this work could prove useful in this particular context. Indeed, when 
studying problems involving the time-evolution of the Hamiltonian by an arbitrary 
variation of the "magnetic field" g(t), it is possible to write exactly the time-evolved 
wavefunction using a dynamical Ansatz [37] 

M 

m)) oc |Ai(*)...A M (*)> = n B ^)) 1°) > ( 42 ) 

i=i 

where a model-dependent set of classical equations of motion is obeyed by Aj(t): 



^ = f*(\ 1 (t)...\ M (t) 1 g(t)). (43) 

For an initial state |Ai(0)...Aa/(0)) which is also representable as | / ui(0)... / Ujv-A/(0)) 
using the alternative pseudo- vacuum one can derive a set of classical equations of motion 
for both representations. It is therefore possible to find, at all times, two representations 
of the true time-evolved wavefunction, i.e. \ip(t)) oc |Ai(£)...Am(£)) \iii(t)...^N-M(t)). 
Since the time-dependent state is no longer writable as a solution to a static Bethe 
equation, Slavnov's determinant would not be available to compute expectation values. 
However equation ( 13"!]) still provides the time evolution of the expectation value of 
observables: 



/ , / A, Q ±,z I , /,m (^1 (*)••• ^iV-Af(*)| Si ,Z \Xi(t)...X M (t)) 

{ms ' m) = <^(i)...^-M(i) m^m (44) 

in terms of simple N by N determinants (or a sum of them for S z ). 

We do not claim here any superiority of the proposed A (e^)- dependent determinants 
over the usual Izergin ones ff22l) . We simply want to draw attention to the fact that, 
in this context, form factors can, in principle, be written as partition functions which 
provide simple formulas valid at any time. 
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7. Conclusions 

In this work we studied Gaudin models realized in terms of spins of finite magnitude 
whose spectrum is bounded from above and below such that the Algebraic Bethe 
Ansatz can be carried out using two distinct quantization axes. We showed that the 
correspondence between both representations of its eigenstates is remarkably simple in 
terms of the set of variables A(e,) directly related to the eigenvalues of the conserved 
operators. We derive a determinant representation of domain wall boundary condition 
partition functions written in terms of the variables A(e$). By mixing the two possible 
representations it was then possible to write overlaps and local spin raising (lowering) 
form factors as such a partition function, making them writable in terms of the proposed 
determinant. Finally, we also point out how these ideas can find direct applications in 
the numerical treatment of certain out-of-equilibrium problems. 
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